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Abstract 



, This paper considerers the problem of computing the value of a solution of the 

heat equation at a given point inside a bounded domain after the initial time. It 
is assumed that the initial value of the solution inside the domain (possibly in a 
part of the domain) is known; the boundary value and the normal derivative on 
a part of the boundary of the domain over a finite time interval are known. Two 
analytical formulae for the problem are given. Both formulae make use of a special 
fundamental solution having a large parameter of the backward heat equation. 
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O 1 Introduction 

^ i In this paper we consider the following inverse problem. Assume that we have a known 
^ I heat conductive body and know the initial temperature inside the body (say constant). 

We are in the situation that we can not access the whole boundary, however, want to 
know the time evolution of the temperature at a given point inside the body (possibly 
close to an inaccessible part of the boundary) after the initial time. How can one know 
it from the temperature and heat flux on an accessible part of the boundary of the body 
for a finite observation time? 

This is a typical inverse and ill-posed problem and appears in many areas of engineering 
and medicine. The aim of this paper is to develop an analytical approach to the problem. 
We consider the simplest, however, important case: the heat conductive body has a 
known isotropic and homogeneous heat conductivity; there is no heat source or sink 
inside the body. In this case, using Fourier's law, after a scaling, one may assume that 
the temperature at a given point and time inside the body satisfies the heat equation. 
Then the problem is formulated as follows. 

Let Q C R"(n = 1, 2, 3) be a bounded domain with a smooth boundary and < T < 
oo. Let u = u{x,t) satisfy 

dtu = Au mnx]0, T[. (1.1) 
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Given non empty open subsets F C dQx]0, T[ and U G (typically U — Q) find a formula 
for calculating M (a;, t) inJlxjO, T[ from the data O)forx e U and {u{x,t), du/ dv{x,t)) 
for (x, t) G Fx ]0, T[. Here v denotes the unit outward normal vector field to dVt. 

We do not assume that u{x,t) for {x,t) G {dQx]0, T[) \ F and u{x,0) for x E Q\U 
are known. These are considered as unknown input or output that one cannot control. In 
this paper we give two solutions to the problem. The consequence of the results stated in 
Sections 4 and 5 is the following. 

Let (xo,to) e nx]0, T[,uj e and c> 0. If x{0}, fix {T} and (^fix ]0, T[)\ 

F are contained in the half space {x t)'^ ■ uj{c) < {xq t^y ■ u{c) — 5 in R"+^ with a 5 > 
and uj{c) = {l/VT+~^){cuj - l f{e S"), then there is an explicit formula to calculate 
u{xo,to) from the data. See Figure 1.1 for the configuration of U, F and the half space 
{x t)^ ■ a;(c) < {xq to 



Here we explain our approach together with the construction of the paper. 
1.1 First approach 

First we observe that the complex exponential function e^'^~*(^'^), z G C" satisfies the 
backward heat equation {dt + A)w = in R""*"^. Choosing a suitable z depending on a 
large parameter r, direction uj and c (see (4.6)), we see that the real part of the phase 
iunction x-z—t{z-z) has the form t\/1 + {x t)^-uj{c). This yields that if (a; t)^-uj{c) > 0, 




u 



Figure 1.1: Configuration. 
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then |e^'^~*(^'^)| — > oo as r — > oo; if {x t)^-cu(c) < 0, then |e^'^-*(^'^)| — > as r — > oo. 
This means that the asymptotic behaviour of e^'^~*(^'^) divides the whole space time into 
two parts whose common boundary is the hyper plane {x t)^ ■ a'(c) = 0. 

Second we choose D C R""*"^ with D G Qx]0, T[ and {xq, tg) G dD in such a way that 
for all p smooth on D and for constants jj, and C independent of p 

lim r'*e-^°-^+*°(^-^) / e'''~*^''^p{x,t)dxdt ^ Cp{xo,to). (1.2) 

r >oo J D 

The examples of the choice of D are given in subsections in Section 4. 

Third we construct a special solution of the backward heat equation with a special 
inhomogeneous term: 

dtv + Av + e'^-'-'^'-'^xoix, = in R'^+^ (1.3) 

The point of the property of v is: the growth order of a suitable norm of e~^'^+*(^'^)v as 
T — > oo is at most algebraic. For the construction of v we basically follow a Fourier 
transform method for a reduced equation to (1.3) which has been used in [13] since it is 
constructive. It starts with introducing a special fundamental solution denote by Gz for 
the reduced equation. This together with necessary estimates are described in Section 3. 

Finally, integrating the equation (1.1) multiplied by v over Qx]0, T[ and applying an 
integration by parts formula described in Section 2, we know that the quantity 

^n^-xo-z+toiz-z) r e^-'-^'^'-'^u{x,t)dxdt 
Jd 

which is coming from the inhomogeneous term in (1.3) multiplied by -y-Me"^"'^"'"*"^^'^', is 
divided into two parts: the first part consists of the data on F and U x {0}; the second 
unknown data on {fl \ U) x {0}, {dflx]0, T[) \ F and Q x {T}. However, from the 
property of v mentioned above and the assumption on the configuration of {Q\U) x {0}, 
[dQx ]0, T[) \ F and Q x {T} relative to the hyper plane {x t)'^ ■ uj{c) = {xq to)'^ ■ uj{c) we 
know that the second part tends to as r — > oo. Combining this with (1.2) for p = m 
on D, one gets a formula to calculate u{xo, to) from the data on F and U x {0} only. See 
Section 4 for the precise description of the formula. This approach can be considered as 
an application of the enclosure method [8] to the heat equation (1.1). 

1.2 Second approach 

In the final section we give an integral representation formula of Gz that is important for 
calculating v. Moreover as a byproduct we present a Carleman type formula for the heat 
equation, which instead of v, makes use of e^'^~^'^^'^^Gz{x.,t) directly which is a special 
fundamental solution of the backward heat equation. This approach is classical and the 
formula can be considered as an extension to the heat equation of Yarmukhamedov's 
formula [15] which has been established for the Laplace equation. 

1.3 Sideways heat equation 

There are extensive mathematical studies in the case when the one dimensional heat 
equation is appropriate. This case is also referred as the sideways heat equation. Carasso 
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[2] considered the heat equation in the half hne 

dfU — Uxx in]0, oo[x ]0, oo[ 

with the initial condition u{x, 0) = for < x < oo. He gave a Tikhonov regularization 
procedure that yields an approximation to u{x, for < t < cxd at an arbitrary fixed 
< X < 1 from noisy u{l^t) for < t < oo. Note that by solving an initial boundary value 
problem for the heat equation in 1 < ,t < oo, one gets the heat flux u.j.{l,t) from u{l^t) 
and thus his problem is reduced to one- dimensional version of our problem. Motivated by 
Carasso's work Levine [11] considered a radially symmetric solution of the heat equation 
in higher dimensions and established a Tikhonov regularization procedure. Their methods 
are based on an explicit integral representation for the solution of the heat equation and 
not time local in the sense: to determine the value of the solution of the heat equation 
at a fixed point in the body one needs the data for all the time (T = oo). In contrast 
to their methods our method does not make use of any representation of the solution of 
the heat equation and needs data only on an appropriate finite time interval depending 
on the location of the point and time where and when we want to know the value of the 
solution and covers fully multi dimensional cases. For possible applications and numerical 
studies of the sideways heat equation, see [1, 4] and references therein. 



2 A weak formulation of the direct problem and in- 
tegration by parts formula 

In this section we describe what we mean by a solution of (1.1) together with an integration 
by parts formula. We follow the formulation of the direct problem described in [3] . 

Given p E L°^{dn), fo E L'^{0, T; (H^Q))') and ho E L^{0, T; H'^/^idQ)) we say that 
u E W{Q, T;H\n), {H^n))') satisfy 

dtu-Au^fo inQx]0, T[, 

(2.1) 

Wu ■ u + pu^ ho ondflx]0, T[ 
in the weak sense if the u satisfies 

< u'{t),v > + / Vu{x,t) ■Vv{x)dx + / u{t)\aa ■ v\da pdS 

=< fo{t),v > + < ho{t), v\9n > in (0, T), 

in the sense of distribution on (0, T) for all v E H^{Q). 

Note that by Theorem 1 on p.473 in [3] we see that every u E W{0, T; H^{n), (H^Q))') 
is almost everywhere equal to a continuous function of [0, T] in L^(f2). Further, we have: 

w(o, T; H\n), (H\n)y) ^ c\[o, t]-, L'{n)), 

the space C°([0, T]; L'^{fi)) being equipped with the norm of uniform convergence. Thus 
one can consider u{0) and u{T) as elements of L^{Q). Then by Theorems 1 and 2 on 
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p. 512 and 513 in [3] we see that given uq G L'^{Q) there exists a unique u such that u 
satisfies (2.1) in the weak sense and satisfies the initial condition u{0) = Uq. 

Given p E L°^{dn), /i E L^O, T; {H^Q))') and h E L'^{0, T; H-^''^{dn)) we say that 
a V e W{Q, T;H\n), {H^n))') satisfy 

dtV + Av = /i inQx]0, T[, 

(2.3) 

Vv ■ 1^ + pv = hi ondQx]0, T[ 
in the weak sense if the v satisfies 

<v'(t),v>- Vu(x,t) ■Vv(x)dx - / u(t)\on ■ v\on pdS 
Jn Jan 

=< fi{t),v > - < hi{t), vlon > in(0, T), 

in the sense of distribution on (0, T) for all v E H^{Q). 

Proposition 2.1. Let u E W{0, T; H^{n), {H^{n)y) satisfy (2.1) in a weak sense; let 
V E W{0, T; H^{fl), {H^{Q)y) satisfy (2.3) in a weak sense. Then the formula 

[ (< hi{t),u{t)\9n > - < ho{t),v{t)\an >)dt 
Jo 

= / u{x,0)v{x,0)dx- I u{x,T)v{x,T)dx+ ^ {< fi{t),u{t) > + < fo{t),v{t) >)dt, 
Jn Jn Jo 

(2.4) 

is valid. 

Proof. Substituting v — v{t) into (2.2), we have 

< u'(t),v(t) > + [ Vu(x,t) ■Vv(x,t)dx 

Jn 

+ f uit)\en-v{t)\9npdS (2-5) 

Jan 

=< fo{t),v{t) > + < ho{t),v{t)\9n > a.e.in(0, T). 
Similarly we have also 

< v'{t),u{t) > - [ Vv{x, t) ■ Vu{x, t)dx 

Jan 

- I v{t)\Qn-u{t)\QupdS (2-6) 
Jan 

fi{t),u{t) >-< h{t),u{t)\9n > a.e.in(0, T). 
Taking the sum of (2.5) and (2.6), we have 

< u'{t),v{t) > + < v'{t),u{t) >=< foit),v{t) > + < hoit),v{t)\9n > 

+ < fi{t),u{t) > - < hi{t),u{t)\dn > ■ 
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Integrating both sides of this equation over the interval (0, T) and using the 
(Theorem 2 on p.477 in [3]) 



/ < u'{t),v{t) >dt+ [ < v'{t),u{t) > dt 
Jo Jo 

— / u{x,T)v{x,T)dx — / u{x,Q)v{x,Q)dx, 



we obtain (2.4). 
□ 

In particular, let u e W{0, T;H\n), (H^n))') satisfy 

dtu - Am = inl^xjO, T[, 

V-u ■ 1/ + fm^ ho ondD,x]0, T[ 
in the weak sense and v G W{0, T; H\Q), (H^Q))') satisfy 

dtV + Av^fi inQx]0, T[, 

Vv ■ u + pv = hi ondQx]0, T[ 
in the weak sense. Then (2.4) gives the formula 

/ (< hi{t),u{t)\3n > - < ho{t),v{t)\3n >)dt 
Jo 

= r < fi{t),u{t) > dt+ [ u{x,0)v{x,0)dx- [ u{x,T)v{x,T)dx. 
Jo Jn Jn 

3 A solution of the equation dfV + Av + / = 

In this section given / = f{x,t) we construct a special solution of the equation 

-dtv - A-u = / in R''+^ 

Note that — 9t — A is the formal adjoint for dt — A. 
Given a complex vector z e set 

v{x,t) = e''-'-'^'-'^w{x,t), f{x,t) = e"-"-*("")^(a;,t). 

Then (3.1) becomes 

—dtw — 2z ■ Vw — Aw = g. 
Taking the Fourier transform of both sides, we obtain 

{-ir)-2iz-^+\^\^)w{^,r))^g{C,v)- 
This motivates us to study the meaning of 1/Pz{^,rj) where 
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Let 

z ^ a + ib, a, be R". 

Since 

PM,r])^\^ + b\^-\b\^-i{r^ + 2a-0, 

we have: PzH-, r/) = if and only if |^ + 6| = |6| and r] -\-2a ■ ^ = 0. Therefore the set 

{{^^ff) G R""*"^ \ Pz{S,^'ri) = 0} is a compact set and forms a submanifold of R""*"^ with 
codimension 2 provided 6 7^ 0. This yields that for z £ with Im^ 7^ the function 
1/Pz{i,ri) defines a tempered distribution on R'^^^. 

Definition 3.1. Define Gz{x,t) as the inverse Fourier transform of 1/Pz{^,r]): 

In this section we study the convolution operator g * g acting on the set of all 

rapidly decreasing functions on R""*"^. 

First we restrict ourself to the case when Re^; = 0. More precisely, given u; e S'^"^ 
and T > set 

Z — ITUJ. 

We set F-,-{x, t) — Gz{x, t) for this z, that is 

F (r f)= ^ f r'^^O+itr, ^^^'n 

' ' (27r)"+i J -irj + 2Tu;-^ + \^\^' 

Let us study the property of Ft-{x, t). For the purpose we employ the argument done in 
[7]. The points are: a relationship between the operators F^- * ■ and Fi * ■ ; an estimation 
of a scaling effect on weighted L^-norms; an weighted L^-estimate for the operator Fi * • . 

3.1 Scaling laws 

It is easy to see that Ft-{x, t) satisfies the scaling law: 

VA > Fxr{x, t) = X^'FriXx, XH). (3.4) 
Given a distribution g{x,t) define 

gx{x,t) ^ g{Xx,XH), A > 0. 
It is easy to see that {gx-i)x — g. From (3.4) we have 

Ft * g{x, t) = j Ft{x -y,t- s)g{y, s)dyds 
= r" I Fi{t{x - y),T\t - s))g{y, s)dyds. 

This yields 

{Fr*g){T-^x,T-H) =t'' j Fi{x - Ty,t - T^s)g{y, s)dyds 



^ J Fi{x-y,t-s)g{T ^y,T '^s)dyd& 
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That is ^ 

{Fr * g)r-^ = :^Fi * (g^-i) 

or equivalently 

Fr*g = ^{Fr*{g,-i)}r. (3.5) 

This yields also 

D:D^F, * = r-2+l"l+2'^ {/^^A^'i^i * (^r-i)}, • (3.6) 

Given s G R we denote by L^(R""'"^) the set of all tempered distributions g — g{x,t) 
that satisfies (1 + |xp + t'^y/^g e L^{W+^) and set 



\g{x,t)f{l + \xf + tydxdt^ 



1/2 

, s e R. 



Note that the set of all rapidly decreasing functions on R""^-*^ is dense in L^(R"'^^). 
Given R > let r > R. Set 



Then we have 



C{R) ^mm{R\R\l}{> 0). 

T 



Since 



= j \g{TX,TH)Wl + \x\'^ + tydxdt 

If ( \x? i^y 



from (3.7) one gets 

C{RY'^ 



M^<75iT(;^NU ^<o. (3.8) 



For T > Rwe have 

1 I 2| |2 I 4j.2 4/ -4 I -2| |2 , j.2\ 

1 + T \x\ + T t — T [T +T \x\ +t ) 



^''(i' + ^-W^^) (3.9) 



Let s' > 0. Since 
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from (3.9) we obtain 

^2s'+(n+2)/2 

\\9r-4s' < -c(RyJ^^\3\\s'. (3.10) 
3.2 Weighted L^-estimates 

Lemma 3.1. Let — 1 < 5 < 0. Given a rapidly decreasing function g on R""^-*^ the 
tempered distribution Fi * g belongs to LKR""*"^) and there exists a positive constant Cs 
independent of g and cu such that 

WD^D^F, * g\U < Cs\\g\\i+s, \a\ + 2/5 < 2. (3.11) 
Proof. For z — iu, we have PzH-, rf) — \^ -\- — 1 — ir}. Let |^| > 8. We have 

= ler + 4(e-(^)|er + 4(e-a;)2 + ry^ 



Next let \i\ < 8 and |(C, r/)| > SVTTW. We have 

|^|2 > 82(1 + 8') - lei' > 8^ > lel^ 
This yields \r]\ > (l^p + \r]\)/2 and thus one gets 

\m,v)\'>v'>Q\\cn\v\r. 

Therefore it holds that, for all {^,r]) e with |(^, 77)] > SVTTW 

\m,v)\>l{\e+\v\)- 

Using this inequality, a local representation of 1/Pz{^,f]) in each neighbourhood of some 

zero points of Pz{i,rj) and Lemma 3.1 in [13], we have the desired conclusion. 

□ 

A combination of (3.5), (3.6), (3.8), (3.10) and (3.11) yields 

Proposition 3.1. Let — 1 < 5 < and R> 0. For all rapidly decreasing functions g on 
R**"*"^ and t > R we have 

WD^D^Gir^ * g\\s < -^n-B^urhh+S: H + 2/3 < 2. (3.12) 
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Proof. First consider the case when \a\ = (3 = 0. We have 

\\Ft * gWs = ^\\{Fl * (^r-0}r||<5 

^ C{Ry/'Cs 

- ^25+2+(n+2)/2 Ml 1+5 

- ^25+2+(n+2)/2 (^(^)(l+5)/2 

Since we have (3.6), a similar argument yields (3.12) for the case when \a\ ^ or (5 ^ 0. 

□ 

For general z the following property is the starting point. 
Proposition 3.2. For all z e C" with Imz ^ Q we have 

G^{x, t) = G^j^^{x - 2tRez, t). (3.13) 
Proof. Let z — a + ib. From (3.3) we have 

^-l^'^J (27r)"+W ' -i(77 + 2a-0 + 26-e+|eP' 

Then change of variables 

77' = 77 + 2a-e, = e 

yields 



^ ' ' 27r y 



(27r)"+i J -irj' - 2i{ib) ■ f + |f | 

Since 

i(a; ■ ^') + it(77' - 2a • ^') = i(a; - 2ta) ■ ^' + itr]', 

we obtain the desired formula. 

□ 

Now given a real vector c set 

gc{x,t) = g{x - tc,t). 
From (3.13) we see that, ior z — a + ib 

{G^ * g){x, t) = J G^{x -y,t- s)g{y, s)dyds 
^ / ^ib^^ - y - 2(i - s)a, t - s)g{y, s)dyds. 
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This yields that 

{Gz *g){x + 2ta, t) ^ j G.jj{x -y + 2sa, t - s)g{y, s)dyds 

= j G^fj{x -y',t- s)g{y' + 2sa, s)dy'ds 

and thus we have 

{G,*g).2a = G.ij*{g.2a). (3.14) 
Since {gc)-c = g, from (3.14) we obtain 

Gz*g^ {G^lJ * {g-2a)}2a (3.15) 

and also 

D^(G,*y) = (D,"G.5*y_2a)2a, 



" / (9 \ 
Dt{Gz *g) = - -^G-Ij * g-2a ^aj + {DtG.jj * g-2a)2a- 



(3.16) 



Remark 3.1. The equation (3.15) corresponds to the simple fact: a function w{x,t) 
satisfies the equation (3.2) if and only if the function w{x, t) = w{x + 2ta, t) satisfies the 
equation 

—dtw — 2ib ■ Vw — Aw — g{x + 2ta, t). 

Now we give an estimate for * g- 

Theorem 3.1. Let — 1 < 5 < and R > 0. Let z — a + ib. For all rapidly decresing 
functions g on R'*+-^ and b ^ with \b\ > R we have 

\\D:Gz * ^11^ < C{R,5) (^l + |aP + |a|) |6|l"l||^|U+i, \a\ < 2 

\\DtGz*g\\s < C{R,S) (^1 + + kl) (2|a||&| + I^Dlbllm- 

Proof. It suffices to consider only the case when a ^ 0. Let \a\ — 0. Set / = G^i5*(fl'— 2a)- 
From (3.15) we have 

l|G'.*^||Hll/2all' 
= J \f(x-2ta,t)\^{l + \x\'^ + tydxdt 

= / \f{y,tmi + \y + 2ta\' + tydydt. 



One can write 



fy 

\y + 2taf + t^ ^ A 

V t 



(y 
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where 

I In 2a \ 

A = 

\ 2aF 1 + 4|a|2 / 

It is easy to see that the eigenvalues A of yl coincides with the roots of the equation 

A2-2(l + 2|anA+l = 0. 

Solving this equation, we have 

A = (1 + 2\af) ± 2|a|^l+ |a|2. 
Since the minimum eigenvalue has the form 

1 



(l + 2|ar) -2|a|^l + |a 
one has 



(Vl + |aP + |a|f' 



: 




: 









(Vl + |aP + |a|f 
Since — 1 < 5 < 0, from (3.17) we have 



<^ ■ \<U\ + \a\^ + \a\)\\yf + t\ (3.17) 



-2<5 



ll/2all6<(Vl + l«l' + l«lj 11/11^ 

and 

2(1+5) 



2 / / /i , I„I2 , l„l \ Il„ll2 

5+1- 



These together with (3.12) yield 



< C{R,S) (Vl + |aP + |a|)"'' |b_2all^. 



/ ; N -25+2(1+5) 

<C(i?,5) (Vl + |aP + |a|j ll^ll^^,. 

Other cases also can be proved as above since we have (3.16). 

□ 

Therefore the map g ^-^ Gz * g & LKR""*"^) can be uniquely extended as a bounded linear 
operator of L^_,_j^(R"+^) into L^(R"+^). We denote it by the same symbol. Then we see 
that, given z with Im^; 7^ and g e L^_,_i(R"+^) the the function 

v{x, t) = e^-^-*("-")(G'^ * g){x, t) (3.18) 
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satisfies the backward heat equation with a source term in the sense of distribution: 

dtv + Av + e'-^-^^^-^^g = inR"+\ (3.19) 

Note that e~^'^'*"*^^'^)i'(x, e L^(R"''"''^) and this v is unique. This is a consequence of 
Theorem 7.1.27 of [6] and the facts that the set of all zero points of -P(^, ^7) has codimension 
2 in R"^"*"^ and — 1 < 5 < 0. See also Corollary 3.4 in [13] for this type of argument. 
Remark 3.2. Hsich[5] developed the scattering theory associated with the operator dt — A 
in R^+^. For the purpose he studied the operator L^(R^+^) 3 f ^ {^/Qz{i,v))'l * f ^ 
L^(R^+^) where q the Fourier transform of a function q on R^+^, z E C"^ and QziC, v) the 
symbol of the operator e~^'^~*(^'^)(9t — A)e^'^"'"*^^'^). Note that e^'^+*(^'^) satisfies the heat 
equation not the backward heat equation. In the paper there is no result related with 
Theorem 3.1. 



4 A computation formula of u in R^^^ with n = 1, 2, 3 

In this section we assume that ho in (2.7) satisfies ho G i^^(0, T]L'^[dVL)) not just ho G 

L2(0, T;H-^i\dn)). 

Let u e W{0, T;H^{n), {H^{n)y) satisfy (2.7) in the weak sense. Let D a bounded 
open subset of R"+^ with D C fix ]0, T[. We denote by xd the characteristic function of 
D. Let V be the distribution given by (3.18) with g = xd- 

Now choose a sequence gj G C^{W~^^) in such a way that gj — > xd in -^5+1 (R"^^) 
as j — > 00. Define 

v^{x,t)^e^-^-'^^-^\G,*g^){x,t). (4.1) 

A combination of Theorem 3.1 and the Sobolev imbedding theorem in R""*"^ we see 
that Vj e C~(R"+^) and {vj|nx]o,T[} is Cauchy in //2,i(Qx]0, T[) = L\0, T;H'^{n)) f] 
H^{0, T;L^{n)){see pages 6-7 in [12]). Since t;,|Qx]o,T[ v\nx]o,T[ in L^{nx]Q, T[), we 
conclude that v e H^'^{flx]0, T[) and this v satisfies 

||L':A^(e-^+*(-%(x,t))|U.(nx]o,T[) = 0{\z\'), \a\ + 2(3 < 2. (4.2) 

By the trace theorem (Theorem 2.1 on p. 9 in [12]) we have Vj\QQx]o t[ — ^ v\dnx]o t[ 
in H'/''^/\dQx]0, T[) = L\0, T; H'/\dn)) D H^/\0, T-L\dQ))- dvj/dv\au^^o,T[ 
dv/du\onx]o.T[ in H'/^-'/\dnx]0. T[) = L\0, T; H'^^dQ)) f] H'/\0, T; L''{dn)); 
Vj{x, 0) — > v{x, 0) and Vj(x, T) — > v{x, T) in H^{VL). 

Note that Vj G 1^(0, T; H^{fl), {H^{d))') and satisfies (2.8) in the weak sense with /i = 
—e^'^~^^^'''^gj{x,t) and hi = dvj/di' + fWj on dflx]0, T[. 
Thus (2.9) yields 

J I + pvj^ u{t)\9n - ho(t)vj(t)\9n^ dSdt 
I I e'''^~^^'''^^gj{x)u{x,t)dxdt+ j u{x,0)vj{x,0)dx - j u{x,T)vj{x,T)dx. 
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Taking the limit j — > oo, we obtain 



(fdv \ ] 

I { —(x,t) + p(x)v(x, t) u(x, t) — ho(x, t)v(x, t) > dSdt 
anx]o,T[ [yai/ / J 

+ / u{x,T)v{x,T)dx — / u{x,Q)v{x,Q)dx 

^_ f e'-'-^'-'Kix^mxdt 
Jd 



(4.3) 



Note that we made use of the fact that: every e -^^(0, T; L^(9Q)) can be identified with 

^{x,t) = (f){t){x) G L2(9(^x]0, T[). 

Divide aOx]0, r[= V U {(aOx]0, T[) \ r} and O = t/ U (O \ [/). Define 



dSdt 



From (4.3) we have 



— / v{x, 0)u{x, 0)dx. 
J u 

Ur) = - [ e'-'-^'-'^^uix, t)dxdt + R 

JD 



where 



R—— I vix,T)uix,T)dx + I vix,Q)uix,Q)dx 
Jci Jn\u 



( dv 



, T^ix^ t) + pix)vix^ t) u(x, t) — ho(x, t)v(x, t) > dSdt. 

9f^x]0,T[\r \ \dl' V / V / y V / V / V /I 



Given c > and u) E S"" define 



(4.4) 



(4.5) 



uj{c) = 



CUJ 



\ -1 



and for r with c^r > 1 set 



z ^ < 



CT [uj + —cu^ ] , cu^ e S 



n-l 



if n = 2, 3 



cr I 1 + iW 1 — ^ uj, uj — {1, —1} if n = 1. 



(4.6) 



One can write 



Re{x-z- t{z ■ z)} = rVl + 



x 



uj{c) 
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since Rez ■ z — t. 

In our method the concept introduced in the following plays an important role. 
Definition 4.1. We say that D is visible at {xo,to) E R""*"^ as r — > oo from the 
complex direction z given by (4.6) if there exist /i > and constant C ^ such that for 
allpe C°°(D) 



r >oo y £) 



(4.7) 



The constant C is unique if it exists. 

Theorem 4.1. Let {xo,to) G r2x]0, T[ be an arbitrary fixed point. Assume that T > 0, 
Lo, r and U satisfy the following conditions: 

I X 



sup 



sup 



I X 

V 



a; c < 




■uj{c); 



uj{c); 





( ] 




/ 




sup 




■uj{c) < 






{x,t)€{anx]o,T[)\r 






\ 


to ) 



u;{c). 



(4.8) 



(4.9) 



(4.10) 



Assume that D with D C r2x]0, T[ is visible at {xq, to) from the complex direction z given 
by (4-6)- Let v be given by (3.18) with g = xd- Then we have 



u 



(xo,to) = -4 lim 7-'^e-^o-^+*°(^-^)7(T), 

(J T >00 



(4.11) 



where 
Kt 



dv 



^ ^x, t) + p{x)v{x, t)^ u{x, t) — ho{x, t)v{x, t) I dSdt — J v{x, 0)u{x, 0)dx. 



Proof. Assumptions (4.8), (4.9), (4.10) together with the trace theorem, (4.2) and (4.5) 
ensure that |-7-Me~^o'^+*o(^'^)i?| is decaying as r — > oo. Since the heat operator 9^ — A is 
hypoelliptic, we know u e C°°(!D). This together with (4.4) and (4.7) yields (4.11). 

□ 

The formula (4.11) can be considered as an application of an idea in [8, 9] that was 
originally developed for the Cauchy problem for the stationary Schrodinger equation. 

So the problem is reduced to: how to choose D that is visible at {xo,to) from the 
complex direction z. In the following subsections we consider this problem. 



4.1 The case n = 1 

Let 6 > 0. We denote by D{xo,to,uj{c), 6) the inside of the triangle with vertices P — 
{xo,to), Pq — {xq — {S/c)y/l + c^u, to) and Pi — {xo,to + + ) in the space time 
R^+^. The two points Pi and P2 are located on the line (x, t)'^ • u;{c) — (xq, ^o)^ ' <^(c) — 
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In [10] we have already known the following. For the proof see the proof of theorem 
2.3 of in [10]. 

Proposition 4.1. Let D — D{xo,to,oo{c),5). If p e C^{D), then 
lim 2(cT)Ve-^°-^+*°(^-^) / e^-'-'^'-'^p{x,t)dxdt^ t\ Pi - P ol^pjP) 

Therefore this D is visible at {xo,to) from complex direction z. The constants p and C 
in (4.7) are given by /x = 3 and 

(4.12) 

4.2 The cases when n = 2, 3 

The cases when n = 2, 3 start with describing the following which is easily derived by the 
proof of Theorem 2.2 and Lemma 4.1 in [10]. 

Proposition 4.2. Let n >2. Let D C R""*"^ be a finite cone with a vertex at P = {xo,to) 
and a bottom face Q ^ (/) that is a bounded open subset of n-dimensional hyper plane 
{x tf ■ u{c) = {xq tof ■ uj{c) -5. Ifpe C°'%D) with 0<9<1, then 

lim l(cr)"+ie-^°-^+*(^-^) / e'^ '-^^' '^ p(x,t)dxdt ^ Kdp(P) 
T — >oo n\ JD 

where 



Therefore if Kd 7^ 0, then constants p and C in (4.7) are given hj p = n-\- 1 and 

C-"^,. (4.13) 

However, it is not easy to show that Kd 7^ for D with general Q. In the following we 
specify Q and show that Kd ^ 0. 



4.2.1 The case when n = 2 

Let 6 > 0. Choose arbitrary two points ,Ti, X2 on the line x-u = xq-uj—{6/c)\^1 + in such 
a way that the orientation of the two vectors u, Xi —X2 coincides with that of the standard 
basis 61, 62. We denote by D{xo, Xi, X2, u!{c), 5) the inside of the tetrahedron in R^+^ with 
the vertices (xo, to), (xi, to), (2^2, ^o) and (xq, + 5VT+~^). We see that the three points 
(xi, to), {x2, to) and {xo, to + 5\/l+?) are on the plane {x t)'^ • a;(c) = (xq, to) ■ uj{c) — 6. 
Therefore D{xo,Xi,X2,uj{c),6) coincides with the finite cone with a vertex at {xo,to) 
and a bottom face Q that is the triangle in R^+^ with the vertices {xi,to), {x2,to) and 
{xo,to + 6VTT^). 
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Prom (4.2) in [10] we obtain the formula 



Kn = Kd^ ■ (0 - 1)"^ 

;it{Kx^.,_0-^)}{(i^,+ix^,)-^)} ^'^ ^ 

where D = D(a;o, 2:2, c<j(c), 5), Ui = v^, U2-, = J^o S'-re the unit outward normal 
vector to the faces of D that are triangles Ai with the vertices {xo,tQ), {xi,tQ), {xQ,tQ + 
6y/l + c^), A2 with the vertices {xo,to): {x2:to), {xo,tQ + 5vl + c?), A3 with the vertices 
{xQ,to), {xi,to), {x2,to); 1!} = {c{LO + iuj-^) - 1)^. 

Since Corollary 4.1 in [10] ensures this Kd 7^ 0, we conclude that D is visible at 
{xq, to) from the complex direction z. Note that 1^3 = (0 — 1)^ and i/i ■ — 1/2 ■ 1^3 — 0. 
Therefore from (4.13) and (4.14) we have the simpler expression 

1(1/3 X t/2) X (i/i X 1/3)1 

{(l^3XI/2)-^)}{(l^lXl/3)-^)}' ^ ' 



Now choose Xi, X2, 6 in such a way that D{xo, Xi, X2, u}{c), 6) C f2x]0, T[. Then we 
obtain the formula (4.11) for D — D{xo, xi, X2, u!{c) , S) , n — 3 and C given by (4.15). 



4.2.2 The case when n — 3 

Let 5 > 0. Choose arbitrary three points Xi, X2 and x^ on the plane x ■ cu — Xq • cu — 
{S/c)y/l + (? in such a way that the orientation of the three vectors a;, x\ — X2, xs — X2 
coincides with that of the standard basis ei, 62, 63. The four points xq, Xi, X2 and 2:3 
form a tetrahedron A in R^. We denote by v the unit outward normal vector field to dA. 
dA consists of four triangles: Ti with the vertices xq, xi and X2; T2 with the vertices xo, 
X3 and X2; Ts with the vertices xq, xi and x^; T4 with the vertices xi, X2 and x^. Since 
u takes a constant vector on each Tj, we denote the vector by Uj. In particular, we have 
1/4 = —CO. 

We denote by D{xo,Xi,X2,Xs,uj{c),6) the inside of the finite cone in R^+^ with a 
vertex {xo,to + 6\/l + c^) and the bottom that is the inside of the tetrahedron in the 
space t — to with vertices {xo,to), {xi,to), {x2,to) and (0:3, to)- Then the boundary of 
D{xo, xi, X2, Xs, a;(c), 5) consists of five tetrahedrons: Ai with the vertices {xq, to), {xi, to), 
{x2 ,to) an d {xo,to + 5\/l + c^); A2 with the vertices (.xo.to), (a^2,to), (.T3 to) an d (xo,to + 
5VT+?); A3 with the vertices (a;o,to), (a;3,to), (a;i,to) and (a;o,to + 5vT+?); A4 with 
the vertices (a:o, to), {x i,to), {x2,to) and (0:3, to); Q with vertices (a;i,to), (2:2, to), (2:3, to) 
and {xo, to + 5^/l + c^). 

We see that D — D{xo, Xi,X2, x^, uj{c), 5) coincides with the finite cone with a vertex 
at {xo, to) and the bottom Q. Let a be an arbitrary constant vector in C^"*"^. Since 

V(,;,t) • (e^-^-*(^-^)a) = {z - t)^ ■ ae^-^-*(^-^), 
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we have 



3 

= V a • {uj 0)^ / e'-'-'^'-'^dSix, + a • (0 - if I e^-'''^'-'^dS{x, t) 

-a-uj{c) I e''-'-'^'-'US{x,t), 

JQ 

where we made use of the fact that the unit outward normal vector to dD takes {i^j Of 
on Aj for each j — 1, 2, 3; —64 = (0 — 1)-^ on A4; — (^(c) on Q. Since a is arbitrary, one 



obtains 



-T ) Jj^e'"'~'^'"^dxdt 



= E ( 0' ) ^ e^-'-'^'-'US{x,t) +(^\^J^ e^-'-'^'-'^dS{x,t) (4.16) 

-uj{c) j^e'-'-^^'-'^dSix^t). 

Since Q is included in the hyper plane {x tf ■ uj{c) — {xq t^f ■ uj{c) — 5, we have 

g-xo-.+to(.-.)^(^) /• e''-'-'^'-'US(x,t) = 0(e-"^^^). (4.17) 
Jq 

We compute the integral 

W)- [ e^-^-'^^-^^dSix,t), j = l,2,3,4. 
On Aj for each j one can write 

where {a, (3, 7) e Aq = {(a, /3, 7) |q! + /? + 7 < 1, a, /9, 7 > 0} and suitable 
linearly independent vectors in R^+^ and satisfy the condition 

ttj ■ u>{c) < 0, bj ■ uj{c) < 0, Cj ■ uj{c) < 0. (4-18) 

Writing Aj = {aj bj Cj) which is a 4 x 3-matrix and tAq = {{a, /3, 7) jo; + /3 + 7 < 
T, a, 7 > 0}, we have 

^-xo-z+to{z-z)j,^^^ = Jdet AjAj I e("<*^+^^^+T^^>(^ '^^^ dad(3d^ 

J Ao 

1 [ g(aa,+/36,+7C,)-(c(c.+.v/Ml7^c.^) -iVdadpdj. 



JtAo 
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Thus together with (4.18) yields 

(a,-t?)(6,-t?)(c,-t?)' 

where ■d — (c(a; + iuj-^) — 1)^. 

Prom (4.16), (4.17), (4.19) we obtain 

hm ( ^ 1 T-3e-^o-^+*o(^-^) / e^'-'-^^^-^Uxdt 
T — >oo \ — r / Jd 



This yields the formula 



lim e-^o-^+*o(^-^V^ / e'-^-^^^-^^dxdt 



Id (a4-i?)(64-t?)(c4-'J?)' 

By choosing p = 1 in Proposition 4.2, one concludes 



2 4 v^lotATAi 



3! (a4-i?)(64-i?)(c4-t?) 



(4.19) 



Therefore D is visible at (xq, to) from the complex direction z and (4.7) is valid for = 4 
and C given by the formula 



,/dct 

(a4-t?)(64-'«?)(c4-t?)' ^ ■ ' 

Therefore (4.11) is vahd for this C and ji — A provided 5 is chosen in such a way that 
D{xq, xi, X2, X3, oj{c), 5) C Qx ]0, T[. 



5 An integral representation of Gz and a byproduct 

It is quite important for us to compute v given by (3.18) with g = xd- In this section we 
give an integral representation of the distribution Gz{x,t) together with 

K,{x, t) = e'^-'-^^'-'^Cix, t), z^a + ib,by^O (5.1) 

which is a solution of the equation dtV+Av+5{x, t) = in R"+^. Using the representation 
of K^, we show that in the hyper space {x t)'^ • (^(c) < with z given by (4.6) is 
exponentially decaying as r — > 00. As a byproduct of this fact we see that yields a 
Carleman type formula for the heat equation. 
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5.1 Representation of Gz 
Proposition 5.1. It holds that 

G,(x,t) 
^ - 2ta) b - \b\H^ 



' (5.2) 

Proof. First we give a representation for Fi{x,t) = Gii^{x,t) where lu = b/\b\. The 
starting point is the following formulae for the Fourier transform: 

Let H{t) denote the Heaviside function. Let Rec 7^ 0. It is easy to see that: 
if Re c > 0, then 

1 



—iri — c 

if Re c < 0, then 

1 



—IT] — C 

Thus Fi{x,t) becomes 



- J H{t)e-^^e-'^'^dt; 
J H{-t)e-'*e-'*''dt. 



FAx t) = - /e^(=^-^)+^*^ 

' ^ (27r)("+i)y -i77-(l-|e + u;P) 

= — V{/ e'^<{-H(t)e^\^+^\'-^^nd^+ [ e'^<{H(-t)e^\^+^\'-^Hdn 
(27r)"- J\i+u>\<i J\^+uj\>i 



(27r)" J\(+u>\<i (27r)" J|^+c^|>i 

This together with a change of variables yields 

Fi{x,t) 

H{t)e-^-- f , H{-t)e 



(27r)" J|4|<i " {2tt)^ J\^\>i 

and thus we obtain 

ei^-^+tFi{x,t) 
Hit) f .... , H{-t) 



J\£\<1 (27r)'' J\£\>1 



(5.3) 



(27r)'^y|^|<i (2%)^ J\£\>i 

Since 



e*"-«e-l«l'X= (\/^) e-^,a>0, (5.4) 
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n „ 

.|2 



we have, for alH < 

V2y^y (27r)"y|^|<i 
This together with (5.3) yields 



(27r)" i|€|<i \2sjT:\t\j 

2y^y 

It follows from (3.3) that 

G^lj{x,t) = \h\'Fi{\h\x,\h\H). 
A combination of this and (5.5) gives 



1^ 
e 4* 



g-ix • b - \b\'^t 



From this we immediately obtain (5.2) since we have (3.13). 

□ 

5.2 Representation of 

From (5.1), (5.2) and the equation 

1 IxP 
X ■ a — tlap + —\x — 2tap = — — , 

it follows that 

e'^ +w,{x,t) 



where 

w,{x, t) = -e^ • « - ^l«r f MV / e^|6|(x - 2ia) • ^ el^'I'l^l V- 
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Remarks are in order, 
(i) Since the distribution 



n \X\ 



is the fundamental solution of the equation dfV + Av — 0, the Wz{x, t) is an entire solution 
of the backward heat equation. Moreover is a smooth function on the whole space. 

(ii) From (5.8) and a change of variables we know that Wz{x, t) — Wz{x, t). Thus Wz is 
real valued and hence we have 



io,(i,() = -e^'''-*l<' 



Note also that Wz does not depend on the direction of the vector b. 
(iii) Write 

/ 1 \ 1 

\b\'^\^\H + i\b\{x - 2ta) ■ ^ ^ t [\b\^ + i—{x - 2ta) J + - 2ta\'^. 

Combining this with (5.6), we can rewrite (5.8) as 

' (27r)- 4|<|6| 

We study more the expression (5.7). 

(i) The case when t > 0. Since 

xa- t\a\^ + t\b\^ = rVl + ^ j " ^(c), (5.9) 
it follows from (5.7) and (5.8) that 

i^,(a;,t) = -e^v^(^*)"-(^)3V / ^z|6|(a; - 2ta) ■ ^^-(1 - ICDI^I ^ (5.10) 

\2'K ) J\(,\<1 

(ii) The case when i < 0. It follows from (5.4) that 

/ ^\b\{x-2ta)-i^\bmH^, 

^\ - 2ta\^/{At) _ f ^t\b\{x - 2ta) ■ Qh\^\et^^_ 




VV J -^l^^i 
This together with (5.7) and (5.8) yields 



Kz{x,t) = e^vW(-*)"-W fMV / e^l&K^ - 2ta) ■ Qbm' - 1)^^^ (5.II) 

\2TrJ J\^\>i 

Remcirk 5.1. Using the well known formula 

/ e'^<dSiO = (27r)^^/V/>r("-2)/2j(„_2)/2(|r^|r), Vr; e R^ 

J\i\=r 

one can rewrite (5.8), (5.10) and (5.11) as one-dimensional integrals. 
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5.3 Exponential decaying of in the hyper space (x 

and a Carleman type formula for the heat equation 

In this subsection first we show that Kz{x,t) is exponentially decaying as r — > oo if 
{x t)^ ■ Lu{c) < and z is given by (4.6). 

Proposition 5.2. Given 6 > we have, as r — > oo 

sup \K,{x,t)\^0 (e-^^/i+^V") . 

{x tf ■ Lu{c) < -S 

Proof. Let {x, t) satisfy {x t)^ -(^(c) < —S. 

(i) The case when t > 0. Prom (5.10) we have immediately 



\K,{x,t)\<Cne-^^'\br. (5.12) 

(ii) The case when t < 0. Wc divide this case into two subcases: (a) |6p(— t) is large; (b) 
|6p(— t) is not large and can be arbitrary small. 

First consider (a). Given R > let t satisfy \b\H < —R. Prom (5.4) and (5.11) we 
obtain 

\K,{x,t)\ < e-^^ (i!)7ei>i'~''^'"'~'^^^ - ^'^'"'^'i^'":^- ^^-^^^ 

Next consider (b). We employ the expression (5.7) and (5.8). Let t satisfy —R < 
< 0. Using (5.9), we can rewrite (5.8) as 

..{X, t) = if ■ -(^)e-l^l'^ (g) " ^^^^^ e^l^l ' 2^«) " ^ el^^l'l^W^. 

(5.14) 

Since t <0 and -\b\H < R, it follows from (5.14) that 

\w,{x,t)\ < Cne~^^^^\b\''e^. (5.15) 

Since cx ■ u < cx ■ ui — t < — 5vT+~?, we have |x ■ u;| > {5 /c)\/l + and thus 
\x\ > (5/c)\/T+?. Using this together with \b\'^\t\ < R, we obtain 

|^|-ngNV(4t) ^ (|^|2|^|)-n|^|2ng-(|6nx|2)/(4|6ni|) 

<\^,-{ms/cni+c'))/m 
- i?" 

This together with (5.7) and (5.15) yields that 

\K,ix,t)\ < Cr. (^e-^>re^ + ^e-(l^l'(Vc)^(l + c^))/(4i?)^ _ (5_lg) 

A combination of (5.13) and (5.16) gives 

\K.M\ < C„ {e-'^>re« (^ + l) + +M;e-(l''l^(^7c)=(l + .^))/(4fi) j . 

(5.17) 



w 
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Now Proposition 5.2 is a direct consequence of (5.12) and (5.17). 
□ 

Remark 5.2. All the derivatives of Kz{x, t) also have a similar property: for each a e Z" 
and /? e Z+ 

^rVT+^s \d^d^K,{x,t)\ 
{x tf -culc) < -S 

is at most algebraically growing as r — > oo. 

As a corollary of Proposition 5.1 and remark 5.2 we obtain a Carleman type formula. 

Corollary 5.1. Let {xo,to) e f2x]0, T[ be an arbitrary fixed point. Assume thatT > 0, 
Lo, r and U satisfy (4-8), (4-9) and (4-10). Let v{x,t) — Kz{x — xo,t — to) for z given by 
(4-6) and u be a solution of (2.7). Then we have 

u{xo,to)^- lim /(r), (5.18) 

T ►OO 

where 

H'^) ~ / I p{x)v{x, t)^ u{x, t) — ho{x, t)v{x, t)| dSdt — J v{x, Q)u{x, Q)dx. 

Remark 5.3. Yarmukhamedov [15] considered the Cauchy problem for the Laplace 
equation in a three-dimensional bounded domain D that is bounded by the plane 2:3 = 
and by smooth surfaces lying in the half-space x^ > 0. He gave a formula for calculating 
the value of the solution at a given point inside the domain from the Cauchy data on 
the portion in ^3 > of dD. For the purpose he made use of a special fundamental 
solution for the Laplace operator which has been introduced by himself in [14] and is 
parameterized by an entire function E{w), w & C satisfying a suitable growth condition 
when |lm w| — > 00. His fundamental solution takes the form for x — (xi, X2, x^) with 
x' = {xi,X2) 7^ 0: 



^ , , I E{x^ + iJ\x'\'^ + u'^)\ du 

^e{x) = ^ / Im ] , 

27r2 Jo \ X3 + i^\x'\'' + u^ J + 

He chose the special E{w) ~ e'^"' with r > 0. Then ^e{x) has the representation 



u^ 



^ I ^ e^""' ri sin(r^|a;'|2 + M2) / , , „ .\ du 

$E(a;) = -— - / / cos (tJ + ) 

27r2 Jo \ v|x'|2 + ^'2 V V I I 



+ u^ 



From this one sees that the Cauchy data of $e(|/ — a;) on the portion in 7/3 = of dD 
for an arbitrary fixed x & D decays exponentially as r — > 00. This fact corresponds 
to Proposition 5.2 and is an evidence that the formula (5.18) can be considered as an 
extension to the heat equation of his formula. 

It would be interesting to find a hidden 'parameter' for Kz like E for ^e- That is: can 
one find a family of special fundamental solutions for the backward heat operator dt + A 
that contains as a special member? In other words, can one find another fundamental 
solution for the backward heat operator that is decaying in one side of a hyper surface 
not plane? 
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